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VISCOELASTIC FLOWS IN A ROUGH CHANNEL: 
A MULTISCALE ANALYSIS 

LAURENT CHUPIN AND SEBASTIEN MARTIN 


Abstract. In this paper, we consider viscoelastic flows in a rough domain (with 
typical roughness patterns of size e <C 1). We present and rigorously justify an 
asymptotic expansion with respect to £, at any order, based upon the definition 
of elementary problems: Oldroyd-type problems at the global scale defined on a 
smoothened domain and boundary-layer corrector problems. The resulting anal¬ 
ysis guarantees optimality with respect to the truncation error. 


1. Introduction 

Many studies investigate the effect of wall roughness on Newtonian flows. In 1827, 
C. L. Navier j26j was one of the first scientists to note that the roughness could drag 
a fluid. Since then, numerous studies attempted to prove mathematical results in 
this direction, see for instance the works of W. Jager and A. Mikelic [20], Y. Amirat 
and co-authors isnu and more recently the works of D. Bresch and V. Milisic |12) . 
Note that all these works formulate the roughness using a periodic function (whose 
amplitude and period are supposed to be small). In a context of more general 
“roughness” patterns, there exists similar recent results, see muH]. All the previous 
work deal with a Newtonian flow, for which the Stokes or Navier-Stokes equations 
are classically considered. 

Much literature research has been devoted to non-Newtonian fluids, in both math¬ 
ematical aspects and applications. It is well known that numerous biological flu¬ 
ids, blood or physiological secretions like tears or synovial fluids, show these non- 
Newtonian characteristics. In engineering applications people are interested in con¬ 
trolling the flows characteristics to suit various requirements such as maintaining 
the fluid qualities in a wide range of temperatures and stresses. Introduction of 
additives lead to non-Newtonian behavior of the modern lubricants for instance. 
Another application domain is linked to polymers, whose non-Newtonian charac¬ 
teristics appear in a wide range of applications such as the molding or injection 
processes. Some particular classes of non-Newtonian models have often been con¬ 
sidered. This includes the Bingham flow or the quasi-Newtonian fluids (Carreau’s 
law [T31 [22] , the power law or Williamson’s law, in which various stress-velocity re¬ 
lations are chosen [53] or mm for mathematical aspects) and also micropolar ones 
|25) . These models, however, consider the fluid as viscous and elasticity effects are 
neglected. The introduction of such a viscoelastic behavior is primarily described 
by the Weissenberg number, denoted We which can be viewed as a measure of the 
elasticity of the fluid and is related to its characteristic relaxation time. One of the 
laws which seems the most able to describe viscoelastic flows is the Oldroyd model. 
This model is based on a constitutive equation which is an interpolation between 
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purely viscous and purely elastic behaviors, thus introducing a supplementary pa¬ 
rameter r which describes the relative proportion of both behaviors (the solvent to 
solute ratio). Considering the Oldroyd model [2^, the momentum, continuity and 
constitutive equations for an incompressible flow of such a non-Newtonian fluid are, 
respectively. 


p[ dtu + u ■ Vu ) — 77(1 — r)Au -|- Vp — div a = 0, 


div M = 0, 


A ( dt(J -|- u ■ Vcr -|- ga(Vu, a) ] + a — 6Aa = 2rir'B{u) 


( 1 ) 


( 2 ) 


(3) 


In these equations, p, p and A are positive constants which respectively correspond 
to the fluid density, the fluid viscosity and the relaxation time. It is important to 
notice the presence of a term 5Acr, 5 > 0 corresponding to a spatial diffusion of the 
polymeric stresses. Usually (he. for the Oldroyd model) this term is deleted, but 
it can be physically justihed: physical effects that can contribute to the diffusive 
process include hydrodynamic interactions [18], particle diffusion [7] and semiflex¬ 
ibility of polymer blends [23]. Thus, the diffusive Oldroyd models have been the 
subject of intense studies related to the understanding of shear-banded flows or 
phase coexistence, see |2S1 ESI ESI EH ED] and, as a consequence, the investigation of 
the mathematical properties of the diffusive model has gained an increasing interest 
recently [HI El EZl 116] . 

Equations ([I])-([3]) make up a system of 10 scalar equations with 10 unknowns; the 
lubricant velocity vector u = (uj)i<j< 3 , the pressure p and the extra-stress symmetric 
tensor a = (<Ljj)i<jj< 3 . The bilinear application pa, —1 < a < 1, is dehned by 

QaiVu, cr) = a ■ W(-u) — W(-u) • a — a(a ■ D(u) -|- D(-u) ■ a) 

where D(m) and W(u) are respectively the symmetric and skew-symmetric parts of 
the velocity gradient Vu. Usually, D(m) is called the rate of strain tensor and W(u) is 
called the vorticity tensor. Notice that the parameter a is considered to interpolate 
between upper convected (a = 1) and lower convective derivatives (a = —1), the 
case a = 0 being the corotational case IS]. Note that taking r = 1 allows us to 
recover various forms of the generalized Maxwell model. Conversely, a Newtonian 
flow is described by choosing r = 0. 

In this paper, we focus on viscoelastic flows in a rough domain (with typical 
roughness patterns of size e 1). We present and rigorously justify an asymptotic 
expansion with respect to £. The development is done at any order, so that we are 
guaranteed to be optimal with respect to the truncation error. We also highlight 
the particular effects of roughness. 

Several relevant questions are not addressed in this article. First, recent works on 
random roughness, see PE9], could make us think that our results can be extended 
to more general cases of roughness. In fact, the construction of our development 
strongly depends on the behavior of solutions of the Stokes equation on a half¬ 
space, whose lower boundary is periodic. The behavior of such solutions must be 
sufficiently decreasing at inhnity to justify our development. Unfortunately, it seems 
that this decrease is only logarithmic in the case of a random boundary (while it 
is exponential in our periodic case). Second, another task related to the regularity 
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of the roughness patterns is not addressed in this paper; what is the behavior of 
the solution when the patterns are not Lipschitz continuous? In particular, what is 
the influence of roughness jump discontinuities over the flow? Finally, the choice to 
make appear a spatial diffusion (5 > 0) in the Oldroyd model could by argue. From 
a mathematical point of view, it is clearly an advantage since we know that such 
diffusion allows to have solutions to the initial problem dl])-®. Nevertheless, even if 
we admit the existence of a smooth solution to the initial problem without diffusion 
(5 = 0), the development proposed with respect to the roughness parameter £ seems 
to be unsuitable. 

The paper is composed of hve sections. In Section [2l we introduce the diffu¬ 
sive Oldroyd model, we precisely describe the roughness geometry and we recall a 
fundamental result; there exists a solution {u,p,a) to the model. In Section [3] we 
introduce the ansatz, which is a formal asymptotic expansion of the solution {u,p, a) 
with respect to the roughness parameter £ <C 1. By identifying the powers of £ in 
such a development, we obtain some elementary problems at any order. Section H] 
is devoted to the mathematical study of these elementary problems; well-posedness 
and properties. Section [S] provides a rigorous justihcation of the asymptotic expan¬ 
sion by analyzing the remainder and deriving error estimates. In Section [6l we show 
that it is possible to effectively determine all the contributions of the ansatz. We 
will notice that regarding each elementary problem and their overlaps, this crucial 
result is not obvious; we prove that the solution of each problem can be built using 
the only previous elementary solutions. 


2. The diffusive Oldroyd system in a rough channel; statement of 

THE problem 

As a matter of fact, the derivation of reduced models is crucial if one aims at 
performing numerical simulations of the flow. However, roughness patterns lead to 
a sharp increase of the computational costs because the mesh of the domain has to 
be built according to the constraints dehned by the roughness patterns. In order 
to avoid such a costly procedure, reduced models can be dehned by considering 
the huid how in the smooth domain (therefore avoiding heavy costs in terms of 
numerical computations by using coarse meshes) and adding a so-called boundary 
layer correction which takes into account the inhuence of the roughness pattern. 
More precisely, for a Newtonian how, we can proceed as follows (see |26i m EQi ng, 
and also dais] for random roughness patterns); 

• at order 0: assume that u is the velocity huid associated to the rough do¬ 
main ujs and Mo is the velocity huid associated to the smooth domain (i.e. the 
domain has been truncated by considering a flat boundary instead of the os¬ 
cillating one and imposing the Dirichlet condition on the smooth boundary), 
see Fig. [1] Then one has 

\\u — uo\\l^{uj^) ^ e and \\u — uo\\m{uj^) ^ \/e. 

• at order 1: in order to counterbalance the error introduced by the truncation 
of the domain, it is possible to introduce a corrector term which leads to the 
justihcation of wall laws. By considering the smooth domain with a slip 
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condition defined as 

= ail + ai2 

= 021 dyU^^'^ + Q!22 dyU^‘^^ 

n® = 0 

on the smooth boundary (the reals aij being coefficients depending on the 
roughness shape), a velocity held u is dehned, leading to a rehnement of the 
approximation. For instance we have (see [20]) 

||m-m||l 2 (^,) < and \\u- u\\hHuj,) < e- 



Figure 1. a) Rough domain, b) Flattened domain: the rough bound¬ 
ary has been truncated by a smooth one. c) Boundary layer domain 
corresponding to a focus on the vicinity of the rough periodic bound¬ 
ary. 


For a non-Newtonian how of the viscoelastic type, we aim at describing the rough¬ 
ness ehects extensively in the following sense: 

• build an asymptotic expansion based upon elementary solutions (i.e. solu¬ 
tions of problems dehned on smooth domains, thus avoiding complex geome¬ 
tries) ; 

• prove in a rigorous way that the asymptotic expansion is valid at any order ; 

• dehne an algorithm associated to an efficient numerical procedure. 

Therefore we consider the Navier-Stokes equations with the Oldroyd model in a 
rough channel. Consider the domain 

■= I (a;, I/) G X M, -eH < 1 / < l|, 

where is the d-dimensional torus, d = 1 or d = 2, and H a smooth periodic 
and positive function. The boundary ujs is denoted and it is composed of two 


connex components: The upper smooth boundary 7 + = x {1} and the upper 
highly oscillating boundary which is denoted 7 ^, see Fig. [2l 
Then we consider the following set of equations: 

Re (m ■ Vm) — (1 — r) Am -|- Vp = div a + f, in Ue, 

divM = 0 , in Ue, 

We (m ■ Vcr-|-Pa(VM, cr))-|-a — D Act = 2rD(M), in 

M = 0 , on 7 ^, 

Bd^a = 0 , on 7 ^. 


( 4 ) 
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Figure 2. Channel with oscillating boundary. 


This system is the non-dimensional version of the system CQi-0. We have intro¬ 
duced the Reynolds number Re, the Weissenberg number We, a relaxation parame¬ 
ter r G [0,1] and the diffusion coefficient D. 

In the next sections, we aim at describing the structure of the solution {u,p) 
with respect to the roughness number £. Before entering into this description, let 
us recall the main mathematical results related to the stationary diffusive Oldroyd 
model. The problem dehned in a strong form can be associated to a variational 
formulation. Then we have (see jl6] for details): 

Theorem 1. Let f G Let Re > 0, We >0, 0<r<l, —l<a<l 

and D > 0. Let us introduce the following constants: 


8|a|C'^,We||/||H-i \/^min(l - r, D) . 

min(l —r, D)2 ’ 4|a|C2^We 




where is a constant which only depends on the domain oOe- 

• Existence. If (7(1) < 1 then, the problem (jl]) admits a variational solution 
(u, a) which satisfies 


2r||V'u||^2 ||(t||^i < Cfiiy 

Moreover there exists a pressure field p G L‘^{Q) such that {u,p) satisfies the 
first equation of the problem (jTj) in the sense of distributions. 

• Uniqueness. The variational formulation of the problem (jlj) admits at most 
one solution if one of the following conditions is satisfied: 

II/I|h-i is small enough; 
b) Re and We are small enough. 

• Regularity. If f is regular then the variational solution is regular and satisfies 
the problem ^ in a classical sense. 

Let us remark that the corotationnal case (namely a = 0) allows us to get rid of 
the smallness assumption on the data. 


3. Asymptotic expansion 
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3.1. Main ideas: ansatz. Let us describe the structure of the solution by using a 
suitable ansatz: 

u{x, y) = uq{x, y) + eUi | j + eui{x, y) + e^U 2 p |) H- 

+ 00 

= X] (uk{x, y) + t/fc (^T, p , 
p{x,y) = po{x,y) + Pi{x,^,^ + ej9i(a;,y) + eP 2 (t, |j H- 

+ 00 

= (pk{x,y) + Pk+i[x,^,^'^ , 

fc =0 

(t{x, y) = (To{x, y) + eSi + eai{x, y) + £^^2 H- 

H-cjo 

= X] {(^k{,x, y) + Sfc p I) ) . 

k=0 

The dehnition of the asymptotic expansion has to be completed by the description 
of the problems satished by the elementary solutions. Then another task consists 
in showing that each elementary problem is well posed. The hnal task consists in 
showing that the level of truncation in the asymptotic expansion is directly related 
to the quality of the approximation of the exact solution. 

In order to identify the elementary problems satished by the elementary solutions, 
we proceed as follows: 

(1) separation of the macroscopic variables {x, y) and the microscopic ones which 
will be denoted (X, Y) = (f, f) ; 

(2) identihcation of terms with the same order with respect to e in the equations ; 

(3) identihcation of terms with the same order with respect to e in the boundary 
conditions. 

Notation 1. ITe will use lowercase letters to denote elements corresponding to the 
real physical domain (x, u, a, u...) and uppercase letters for all that concerns the 
microscopic field (X, U, S, Vl...). 

To complete this subsection, and before getting into the details, let us explain 
the way the asymptotic expansion has been built; at main order, the huid how is 
governed by a classical viscoelastic model with boundary conditions located at the 
hat bottom 1 / = 0 (dehnition of (uq,P o, o'o)). But, of course, the boundary layer has 
been omitted and, in fact, the boundary condition should have been imposed on the 
oscillating boundary instead of the smooth one ; it can be shown that the result¬ 
ing error is of order e (the size of the boundary layer) and, therefore, a so-called 
boundary layer correction is introduced in order to counterbalance the mentioned 
boundary value default (dehnition of (Pi, Pi, Si)). Now if we analyze the approxi¬ 
mate solution dehned as (mq + ^Pi,Po + -Pi, <xo + ^Bi), the equations in the domain 
and the boundary condition at the oscillating boundary are satished by means of 
construction ; unfortunately, the boundary condition on the upper boundary is not 
satished by the approximate solution because of the behavior of the boundary layer 
corrector at inhnity. However the error is of order e and it is is possible to build a 
solution on the smooth domain (dehnition of {ui,pi, Ui)) which counterbalances this 
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boundary value default on the upper boundary. Again, at his step, by considering 
the smooth domain only (in particular, homogeneous boundary conditions are con¬ 
sidered at the lower smooth boundary), the boundary layer has been omitted and, in 
fact, the boundary condition should have been imposed on the oscillating boundary ; 
the resulting error is now of order Thus, if we compare [ueiPeT^e] ~ 
and {ue,P£, Ug) — {uo + eUi + eui,po + Pi+epi, Uo + cSi -t-eui), the error has been de¬ 
creased by an order of magnitude and, besides, the same procedure can be applied by 
introducing a boundary layer correction which counterbalances the boundary value 
default of order at the oscillating boundary. 

In a more general way, the boundary value default introduced at the oscillating 
boundary can be counterbalanced by a boundary layer correction ; the resulting 
boundary value default at the upper boundary can be counterbalanced by a vis¬ 
coelastic flow dehned in the smooth domain. Through this procedure, the resulting 
approximation satishes the equations in the domain, the boundary condition at the 
upper boundary and the error on the boundary condition on the oscillating boundary 
has been decreased by an order of magnitude. 

Let us dehne the two rescaled sub-domains. As a matter of fact, the main flow is 
dehned on the smooth domain whereas, due to the consideration of the roughness 
patterns, the boundary layer is rescaled by the homothetic transformation {X, Y) := 

X y 
e ’ e 

Definition 1. The smooth domain is defined by 

cuo := e T"* X R, 0 < y < l}. 

The upper boundary equals y’*' and the lower boundary corresponds to '-Jq . The 
normal outward unit vector the lower (resp. upper) boundary is n = (0, —1) (resp. 
(0,1 )) on (resp. 

Definition 2. The boundary layer domain is defined by 

= {(x,y) e X R, -H{X)<Y}. 

The boundary {{X,Y) G x R, Y = —H{X)} is denoted T. We denote by 
N := —(VH, 1) the outward vector to the lower boundary. Note that N is not a unit 
vector. 

Notation 2. The usual notation for classical operator of derivation are V, div 
and A. The problems considered in this paper make appear two kinds of functions: 
the first ones, like the velocity u, which only depend on the macroscopic variables 
{x,y), and the others, like the velocity U which depend on {x,X,Y). In the first 
case the classical operators are defined as usual, for instance 

d 

Au = Y^ dl^u + dlu. 

e=i 

In the case of function depending on (x,X,Y), the notations are the following: 

d d 

AU = 9x,U + dfiU and A^U = ^ d^U. 
e=i e=i 
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Notice that the divergence with respect to the variable x is also defined for a function 
U : ^ by 

d 

i=i 

The normal derivative of a function U defined on T is defined by 

d 

d^u = {N ■v)u = 

e=i 

3.2. Elementary problems. Let us define the problem at the main scale: it con¬ 
sists in considering the viscoelastic problem on the smooth domain, i.e. by truncating 
the rough boundary from the initial domain, associated to homogeneous conditions 
on both boundaries. 

■ Main order: 

Re (no ■ Vno) - (1 - r) Auq + Vpo - div((To) = /o + oq, on uq 

div(no) = 0 , on loq 

We(uo-Vao + ga(Vuo,o-o)) + o-Q-DAao = 2rB(uo), on uq 

Uo = 0, on 70" U 7+ 

dndo = 0, on 7 + 

dn(^o = bo, on 70" 

where /o = / is the source term coming from the modelisation, and where the 
constants Oq and bo have to be fixed. 

Remark 1. The choice of constants Oq and bo will be discussed further. Roughly 
speaking, they play the role of a degree of freedom which will be fixed in order to en¬ 
sure the well-posedness of elementary problems and a suitable behavior of elementary 
solutions to be defined. From a practical point of view, we will set ao = bo = 0. 

Imposing the homogeneous Dirichlet condition on the smooth lower boundary 'Jq 
instead of the oscillating one 7 “ is a source of error. Indeed, {uo,Po, Uq) satisfies all 
the equations except the boundary condition at the oscillating boundary 7 “. The 
value at the oscillating boundary can be determined by a Taylor expansion; 

Uo{x,-eH{^)) = ^ d^Uo{x,0). 

k=0 

Since uo{x, 0) = 0, one can check that is of order e, namely 

uo{x,-eH{^)) = -eH{^)dyUo{x,0) + 0{e^), 

and this leading term of order e will be counterbalanced by first the boundary layer 
corrector. 

Remark 2. Note that the other terms (of orders etc.) will be treated and 

counterbalanced in subsequent boundary layer problems. 

Remark 3. The same methodology applies for the Neumann condition related to 
the elastic tensor. 
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As a consequence of the previous remarks, we define the following boundary layer 
problem. 


■ Correction with boundary layer n° 1: 


( 




V 


-(1 - r)At/i + VPi 
div(Pi) 
Pi 

-DASi 


0, 

on 

0, 

on 

HdyUol^-, 

on r 

0, 

on 

(VR ■ Vx)(To|..y- - bo 

on r. 


Remark 4. Let us present the behavior of the solution (the proofs will be given by 
proposition^^ later). 

(1) When Y +cx3, the velocity Ui exponentially decreases towards the constant 
defined as 

Uoo-.= lim [ Ui{;Y). 

y —^■OO 

(2) For bo = 0, the elastic tensor Ei exponentially decreases towards 0. 


At this stage, it can be shown that the approximation satisfies the equations in 
the domain and at the oscillating boundary. At the upper boundary, the Neumann 
boundary condition for the elastic tensor is satisfied because of the exponential decay 
of (To at inhnity. However the homogeneous Dirichlet condition is not satished for 
the velocity held. This is why a so-called main corrector dehned on the smooth 
domain is dehned in order to counterbalance the boundary default introduced by 
the approximation at the upper boundary. 

Dehning the linear operators and by 

:= Re(Mi-VMo -|- Uq ■ Vui), 

:= We(Mi ■ VcTo -F5'a(VMi,(To)), 

/:(c)(cri) := We(Mo ■ VcTi-h 5 'a(Vno,ai)), 

we dehne the following problems: 

■ Main correction n° 1: 


— (1 — r)Aui + Vpi 

= div(cTi) Oi, 

on ojo 

div(Mi) 

= 0, 

on ojo 

-|- C^^\ai) + ai — DA(Ti 

= 2 rD(Mi), 

on uo 

Ml 

= — lim [ Ui , 

on 


Y^oo 

Ml 

= 0, 

on 7o- 

dn(^l 

= 0, 

on y’*' 

dn<Jl 

= bi, 

on 70". 
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■ Correction with boundary layer n° 2: 



-(1 - r)At/2 + VP 2 
div(t/ 2 ) 
U 2 

-DAS 2 

dN^2 


F 2 — cto, 

0 , 

^ 2 , 

-P(VP-V.)d,aoL- 

+(VP-V.)Si- 6 i 


with 


on n 
on n 

on r 

on n 

on r. 


F 2 := div(Si) - Va;Pi + 2(1 - r)Vx • VPi - Re(Mo ■ VPi), 

G 2 '■= 2rD(Pi) + 2D Va; • VSi — We (uq ■ VSi + ga(Vl7i, (Jq)) . 


Remark 5. With a suitable choice of (resp. bi), the problem related to the velocity 
field (resp. constraint field) is well-posed and satisfies the property with exponential 
decay towards a constant (resp. Oj. Moreover, hi does not depend on ai! 


■ Main correction at order k >2: 

f - (1-r)Auk-hVpk = div(ak) + fk + Ok on uq, 


div(uk) = 0 on loq, 

F^^\uk) + + (Tfc - DAo-fc = 2rD(Mfc) + gk on uq, 

Uk = JjUk on 7 +, 

Uk = 0 on 7(7, 

dnUk = 0 on 7 +, 

. dni^k = bk on 7(7. 


where fk and gk only depend on solutions that were dehned previously; 

fc-i 

fk ■ Rg ^ ^ Hi ■ V 'Uk—ii 

i=\ 

k-l 

gk := We'^{ui-V(Tk-i +g{Vui,(Tk-i)). 

i=l 

Higher order correction terms of the asymptotic expansion are dehned as the 
solutions of the following elementary problems. 


■ Correction with boundary layer of order k > 3: 

' -(1 - r)APfc + VPfc = Fk-ak -2 on D, 

div(Pfc) = 0 on D, 

PB(D , Uk = Dk on r, 

—DASfc = Gk on D, 

dN^k = Nk- bk-i on r, 
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where Dk, Gk and Nk only depend on solntions that were defined previonsly. 

Fk := div(Sfc_i) + 2(1 - r)V3, • Vt/fc-i 

+ (1 - r)A^Uk-2 - '^xPk-i + diVj;(Sfc_2) 

k-2 

“Re {{ui + Ui) ■ {S/Uk-i-i + WxUk-2-i) + Ui ■ 'Wuk-i- 2 ) , 


Dk := 

Gk '■= 2rD(17fc_i) + 2r'Dx{Uk-2) ~ ^k -2 + 2D Va; • VSfc_i + D AxTik -2 

/ k-2 


-Wei Uj ■ (VSfc_i_i + Vx^k-i-2) 

V i =0 

+Ui ■ (V(jfc-j-2 + VSfc_i_i + Wx^k-i- 2 ) 

j— 2 ) “ 1 “ I 7 a(VDj+l “ 1 “ xUii (Jk—i —2 F ^k—i—2}^ j 

Nk := {VH-Vx)ak-,L- + {VH-Vx)^k-i 


■ Vx)ak-,-, + dyak-,-,) |,-. 


4. Analysis of the elementary problems: well-posedness and 

PROPERTIES OF THE SOLUTIONS 

The elementary problems related to the bonndary layer correctors with respect 
to the velocity take the following form; 

{ -(l-r)AD + VP = F onD, 

div(D) = 0 onD, 

U = D on T, 

with F e L'^in) and D e L‘^{T^). 

Notation 3. We denote Fj the following Fourier coefficients of a function F defined 
on X (0, +cxd) .• 

F{X,Y) = 

jeJA 

Notation 4. All constants depending only on the domain and on physical constants 
will be considered harmless, there will be denoted G. In the same way, we will use 
the notation Q{Y) to denote any polynomial with coefficients depending only on the 
domain or on physical constant. In particular the quantities G and Q{Y) do not 
depend on the variables x or X. 

Proposition 1. If the following conditions are satisfied 
Fo{Y) <Q(F)e-^, forallY>0, 

Ty{Y) < Q{Y) e-ll^'IR, for all j G \ {0} and F > 0, 
then admits a unique solution {U, P) satisfying VU G L‘^{Q) and P G L‘^{Q). 
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(1) There exists Uoo G such that 
^(r)-t/oo|| + |n(r)| <QiY)e-^, forallY>0, 

+ K\y) < Q{Y) e-^, for all F > 0. 

(2) We have 

\\U{X,Y)-U^\\ + \P{X,Y)\ <Q{Y)e-^, for all X E , Y > 0, 
and, for all i >1, 

\\d^^U{X,Y)\\ + \d^^P{X,Y)\<Q{Y)e-^, for all X E , Y > 0. 

Remark 6. In this Proposition, the vector Uoo can he identified as the limit ofU^iY) 
when Y goes to +cx); 

Uoo ■■= lim [ U E 

Proof. Let us introduce the following decomposition of the vectors U = {U^^\ E 
X M and F = (-Pd), F^)) g R^ x R. Then, we pass to the Fourier transform with 
respect to X. Passing in the Fourier regime, equations satished by {U, P) inside the 
domain Q, in problem can be translated into 



lUIITj^- c/jd + 

— _pd) 

J__^ 

on {Y > 0} 

Vj e Z‘^ 

(5) < 

IldllTf 

— pU) 
j 

on {Y > 0} 

Vj e Z'^ 


ij . cf ' + uf 

= 0 

on {Y > 0} 

Vj G Z'^ 


where , f/® and Pj belong to L^(0,+cx)). Now we solve the Fourier problem 
and describe the behavior of the solution of the Stokes problem. 

(1) Let us discuss the case j = 0. The system reduces to 

=F(id withf/gd) eL\0,+oo), 

-U^^^ with^e L2(0,+oo), 

U^^ =0 with gP^(0,+cx)). 

By integration, this leads us to the following equalities 

C/f V) = 0. 

(7f(K) = C/f(0)-y^ ^ q‘>K)d^jdZ 

^+oo / ^+oo- \ 

+ X U d‘'(0d,)dZ. 

Wiy) = ^"’(0). 

/ + 00 - 

F^^\z)dZ. 




VISCOELASTIC FLOWS IN A ROUGH CHANNEL: 


A MULTISCALE ANALYSIS 13 


By assumption on the source term, we have 




(C)dCdz 


<Q{Y)e 


-Y 


r*+oo 


F^^\Z) dZ 


'Y 


< Q{Y)e 


-Y 


£= 1 , 2 . 


Dehning [/^ := G x R as 

-—~ r+oo / I'+co \ 

f/S’ := 

US’ := (7f(0), 

we obtain, for all y > 0 , 

ih'W) - r/iJ’l + \uP{Y) - c/£'| + ia(r)| < Q(Y)e-'\ 

ld‘’V)| + |(7rV)l + la (r)l < Q{Y)e-''. 


(2) Let US discuss the case j ^ 0. If the source terms and were written 
as then the estimate is proved in |T^ Appendix B], In the present 

case, the source terms only satisfy 

<Q(h')e-lljll'", vy>0, 

\ f ^^\ y )\ < g(y)e-iijii^, vy > o, 

and the proof can be adapted from jT5] by using a comparison principle. In 
particular, we can show that 


UJY) + PAY) <(3(r)e-"'"'', 


for all j G \ {0} and y > 0. 


□ 


The elementary problems related to the boundary layer correctors with respect 
to the elastic constraint take the following form: 

—AE = G on 12, 

Bm'F = N on r. 


pg(ref 2 ) 


with G G F^iVt) and N G 

We note that the problems related to the boundary layer correctors are matricial 
problems: the unknowns have matricial values. The elementary problem that we 
analyze here is a scalar case but can easily be extend component by component. 


Proposition 2 . Assume that 

[ G{X,Y)dXdY = 
J{Y<0} 

If the following conditions are satisfied 

Go{Y)\<Q{Y)e-^, 

^(y) <Q(y)e-idii^, 


ijd 


N{X) dX. 


for all y > 0, 

for all j E 11^ \ { 0 } and y > 0 , 
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then admits a unique solution S G Moreover, we have 

S^(r) <Q{Y)e-^, forallY>0, 

Ej(r) < Q{Y) for all j G \ {0} and Y > 0. 


In particular, we have 

|S(X,F)| < Q{Y)e-^, for all X e T^, F > 0. 

Proof. If the source terms satisfy Go = 0 and, for j > 1, Gj{Y) = then 

the estimate is proved in [TH Lemma 2.2], In the present case, the source terms 
only satisfy 

|^(F)| < Q(F) e-^ and |^(F)| < g(F)e-ll^l"^. 

The estimate on Eq can be obtained by straightforward integration. For j ^ 0, the 
proof can be adapted from H by using a comparison principle. □ 


Applications; analysis of problems Let us recall that the dehnition of the 

boundary layer correction problems PB^^^ need to specify the value of 0^-2 and 6 a:-i- 
Let us hrst describe how to determine 0 ^- 2 . In order to apply Proposition [H we 
need to impose that the source term in the momentum equation of PB*^^\ namely 
F := Fk — afe_ 2 , satishes a sharp decrease for each Fourier mode: 

• Averaging this source term F with respect to X gives 


Fo(y) = / (Ffe(X,F)-afc_2)dX 


jTd 

([ Fk{X,Y)dX- lim [ Fk{X,Y)dX 

\Jjd y->- + cxD Jr^a 




, lim / Fk{X,Y)dX-ak -2 

\Y —>-+00 J '^d 


On one hand, Fk is composed of elementary solutions Pi)i<i<k-i 

which, by induction, satisfy the expected decreasing behavior. We have 



Fk{X,Y) 


lim [ Fk(X,Y)dX 

y "—>-+00 


<QiY)e-^. 


On the other hand, in order to satisfy the assumption needed to apply Propo¬ 
sition [H we impose 


( 6 ) 


at-2= lim / Fi(A',y)dA'. 

Y —>-+oo 


By induction on k, the Fourier coefficients Fj, for j 7 ^ 0, of the source term 
F = Fk- ak -2 satisfy 


Fj{Y) < Q{Y) e-ll^ll’^, Vj e'L^\ {0}, VF > 0. 


Let us now describe how to determine hk-i- In order to apply Proposition m we 
need to impose the compatibility condition between the source term in the Laplace 
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equation of namely G := D and the Neumann boundary term N = 

Nk — bk-i'. 

(7) bk-i= [ Nk{X)dX-D-^ [ Gk{X,Y)dXdY. 

Jid Jy<o 

Besides, by induction on k, the Fourier coefficients Gj of the source term G := D~^Gk 
satisfy 

|^o(V')|| <Q{Y)e-^, VF>0, 

^j{Y) < Q{Y) e-ll^ll^, Vj eZ<^\ {0}, VF > 0. 

5. Error estimates 

5.1. Remainder. The asymptotic expansion truncated at a given order leads us to 
introduce the so-called remainder (IZ, Q,S): 

N 

u{x) = '>^j i^) Y Uj (^x,+ e^Tl{x), 

j=0 ^ 

N 

Pji^) + Pj+i(x,-) + £^Q{x), 

j=0 ^ 

N 

a(x) = (^) Y ~'j Y e^S(x). 

j=0 ^ 

We aim at establishing estimates on the remainder, at any order. 

Applying the Oldroyd operator to the remainder (77., Q,S) and considering the 
properties of the elementary solutions, we get the following set of equations 

• momentum equation, in cUg: 

Re {e^n ■ VIZ + C[Y{Tt)) - (i - r)XTZ + VQ = div(5) + j;, 

• continuity equation, in uje'. 

div(77) = 0, 

• constitutive equation, in ujs'. 

We {IZ-VS + gaiyiZ, S)) + Ci^\lZ) + cP{S)'^ + 5 - DA5 

= 2r'D(lZ) + 

• boundary conditions on the velocity, on 7 “ U q’*'; 

77 = T>W, 

• boundary conditions on the elastic constraint, on 7 “ U q"*"; 
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where operators are defined as follows. 




\k=0 J 

/ N \ 

Y^"i^k + Uk)\ vn, 

\k=0 / 


+e^7l-{yuN + V^U^ 


d^Hn) 


\ fc =0 

) 

\k=0 J ^ 

+e^ IZ ■ (Vctat + VxSat) + Qa i 'VIZ, s^ick + Sfc) j , 

\ k=0 J 


dj^'\s] 


+ e^QaiVuN + V^t/jv, S) 


^ N \ 

Y^'dk + Uk)] vs 

N 1 ^ 

+9a I £^{VUk + VxUk + Vt/fc+l), <5 j 

\fc =0 / 


The source terms are defined by 

= 


Re ■ iVujsi-k + VxUN-k + Vt/AT+i-fc)^ 

-Re Mo ■ (Vmat + VxUn) + (1 - r)AxUN + (1 - r)AxUN 

-VpN - VxPn+1 + div((TAr) + diVa;(SAr), 

Qs = 2r'D(^U]\i') + 2r]D)x{U]\i') — (j^v — SjvTDActjvTD A^-S jv 

—We ^ + Uk) ■ (V(JAr_fc + Vx^N-k + VS^r+i-fc)^ 

-We Mo ■ (Vctat + Va;SAr) 

A-1 / \ 

—We ^ ^ I QaiVuk + VxUk + Vf/fc+i, cjsf-k + T,i\f-k) | 

fc =0 V / 

-We5(a(VMjv + VxUN,(ro), 

and the contributions to the boundary relation are given by 


vd 


vd 


lim [ Uk-Uk\ 
to 

N +00 / 


Y=l/s , 




k=0 i=N—k+l 




N-l 


A/-W = -E 




r^fc+ih=i,u=i/e, 


k=0 


ud = 


N 


+ 00 


.-N 


E^‘ E 


(-eHy 

i! 




k=0 i=N—k-\-l 

Proposition 3. The following estimates hold: 

md<C, WGeldAC. 
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For all i > 0, for all x G we have 

\V',V'*\x)\ < Q (Ij e-i, 

|V'A4<+>(i)| < Q eU, 


v'pi-’WI < Ce'-', 


Proof. 


The estimates for and Qe are obvious. 

By using Proposition [H we know that for each integer k we have 


Uk ( X, - 
e £ 


and, for all £ > 1, 


lim / Uk 

—^OO /qp 


< Q 


e s, for all x G 


V'vO I X, 1 1 


< Q ( - ) e s, for all a; G T 


We immediately deduce that for all £ > 0 and for all a; G we have 


e-l 


We estimate ^ remarking that, using the Taylor formulae, for each 

integer k we can write 


i=N-k+l ■ ^ 

with ^k £ [0, —eH{x/e)]. That implies 

V^-\x,-eH{x/e)) = -g ^ ^ . x, ( t, Cfc)• 

fe=0 ^ + h 

Now, since is a hnite sum, the estimate directly follows from its analysis. 
• Finally, the estimate on is based on the same arguments, noticing 

that the Neumann data, for the remainder, can be reduced to a finite sum 
of boundary terms, namely 

M£^\x, -eH{x/e)) 

^ (-mN-k+i 

= + dyCTk){xXk). 

k=0 ^ '' 

□ 


5.2. Lift procedure. Well-posedness of the set of equations satished by the re¬ 
mainder is obtained by means of construction. Let us point out the fact that the 
remainder satishes a diffusive Oldroyd-type system. 

In this step, we aim at modifying the set of equations by using a lift procedure in 
order to deal with homogeneous boundary conditions and preserve the homogeneous 
incompressibility condition. 
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Definition 3. Let r G C'°°(R) be such that 




0 if 1/ < 0, 
1 


We define T^-bound (in^d iSbound os 


'Abound(a:,I/) = T{y)V'^+\x) + (1 -r(i/))'D( ^(a;), 
•Abound(a;, 2/) = T{y)M^^\x) + (1 - r(2/)) 


and we define IZdiv os a solution of 


f div(7Zdiv) = -div(7?.bound) 

\ '^div = 0 


in Us, 

on 7 “ U 7 '''. 


The lift velocity and constraint fields are thus defined as: 


■ "Abound “1“ '^div: • Abound 

and the lifted velocity and constraint field are thus defined as: 

IZ := IZ — Q '■= Q, S := S — 


Remark 7. The lifted remainder {1Z,S) satisfies a system which is identical to the 
system satisfied by {7Z,S), up to some modifications: 

• the boundary conditions are homogeneous ; 

• the incompressibility condition is (still) homogeneous ; 

• the source terms Qe, ond linear operators have been 

(slightly) modified (see further, paaelW\) 


Remark 8 . The definition of TZdiv is guaranteed by the following result, due to 
Bogovskii [9] (see also cu;-- 

Proposition 4 (Bogovskii). IfTL^ H^{ufij, m>d, is such that 


R = 0, 


then there exists a solution IZ G H^^^{uf) of 

div('R) = H 

IZ = 0 


m Us, 
on 7 “ U 


such that 


\\S/lZ\\H^(uJe) ^ C WTLWh^^^uj^)- 


Thus, the existence of such a lift function IZ relies on the identity TL = 0 with 
TL = — div(7?.bound)- By the Stokes formula, we have 

f div(7?.bound) = [ 11 + f '^ = [ div(7?.) = 0 . 

J^e J 1e J 7 ^ ^ 

Proposition 5. The following estimates hold: 


Proof. 
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• By definition of 7?.bound and iSbound (see El), the estimates 

11 "Abound 11 ^Ce , || inbound || ^Ce 

directly follows from the Proposition [3l 

• Following the Bogovskii inequality (see Proposition 0]) and the Poincare in¬ 
equality, we have 

Il’^div||n2 < C ||7^bound||R2. 

• Since T^uft = 7?.bound + and iSiift = iSbound, the two previous steps imply 
the result of the Proposition O 

□ 

Now estimates for a homogeneous (w.r.t. boundary conditions and incompress¬ 
ibility condition) system have to be established. This is the purpose of the next 
subsection. 

5.3. Estimate on the remainder. The lifted remainder (fR., Q, S) satisfies the 
following system: 

• momentum equation, in Ue'. 

E^Ren- Vn + ('^) - (1 - r)An + VQ = div(5) + 

• constitutive equation, in Us- 

E^ We (n ■VSR ga{vn, 5)) + ci^\n) + cf\s) + 5 - DAS 

= 2rJ}(fR) 

• the homogeneous incompressibility condition, homogeneous Dirichlet condi¬ 
tions for the velocity, homogeneous Neumann conditions for the constraint. 

The linear operators are given by: 

c^in) = ci^\n) + ReE^ + 

= £f)(^)+We£^(^-V5nft + ^a(V'R,5nft)), 

£F(5) = £fH5)+We£^(7?,nft-V5 + ^,(V7?,nft,5)). 

The source terms are given by: 

= .Fg — Re ■ V7?.iift -|- £i'^^(7?.iift) j -1- (1 — r)A7?.iift, 

Qe = 61.-5Hft + DA5uft-We(£P)(7?-uft) + £F(5uft)) 

—E^ We (7?.iift ■ V5iift -|- 5'a(V7?.iift, «Siift)). 

Remark 9. From the PropositionslM andl^we can obtain a bound on the new source 
terms: 

11-^:||l2 < Ce ||^£||l2 < Ce 

Theorem 2. The remainder satisfies: 

\\vn\\L 2 + \\s\\m<CE-\ 
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Proof. Due to the relations between (fR.,S, Q) and (IZyS, Q) 

IZ := IZ — Q '■= Q, S := S — 

using the Proposition [5l it suffices to analyze the error on (TZyS, Q). 

The estimate is then obtained using a classical energy estimate on the system 
satished by (fJZ,S, Q). More precisely, we hrst take the scalar product in L‘^{ue) of 
the momentum equation by 2r'JZ. Next we take the scalar product in of the 

constitutive equation by S. We hnally add the results to obtain 


(9) 


2 r(l 


\\VTir- + 



RHS. 


The term RHS is composed as follow: 


RHS = -2r Ci^\n)-lZ + r R, ■ TZ-Wee^ ga{Vn,S) : S 






- / df\TZ)-S 


CP{S)-S 


<3e-S. 






'LOe 


It is not very difficult to show that the source terms and linear terms of RHS 
can be controlled by the terms on the left-hand side of the estimate (jH])- However 
the quadratic term does not lead in a straightforward way to a suitable estimate. 
Therefore we have to consider a new argument which is based on a hxed point 
procedure. Let us consider the following linearized system, denoted (Rim,): 


e^RelZ VIZ + Cr\TZ ) - (1 - r)A72. + 

= div(5 ) -L j;. 


n+1 


• n -^n+l. 


e^VJeilZ VS +ga{VlZ ,S 


+S 


n+1 


= 2r]D){n^^) + g,, 


and 


n+1 

div7^ = 0. 


~ fl -^ IL 

where [IZ , S ) are given. The idea relies on the following arguments: 


(1) we show that {7Z , «S )„ is bounded in up to smallness assumptions ; 

~ n 

(2) by the Cauchy criterion, we show that the sequence [TZ , S )„ converges in 
H\u,) ; 

(3) we let n tend to -|-oo and show that the limit of {IZ ,S )„ is the solution of 
the system satished by the remainder. The limit still satishes the estimates 
of step 1 . 


Step 1. Using the classical energy estimate, we have 


2 r(l-r) f \\Vlz"^'f+ f ||5"^'|p + D f ||V5"^'||2 

JlU J (jJ J ijJ 


RHS^'^^ 
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The term RHS^"'^ is composed as follow: 


=-2r I ■ 71 +r I 


n+1 


' LOe 


' u)e 


Mi / n -^n+l, ~n+l 

-Wee^ ga{Vn ,S ) : S 




- )-s - cT’is 




' LOe 


'CJe 


+ / ■ «S 


n+1 


'cJe 


We distinguish three types of terms: 
• source terms 


n+1 


J- e + / Ge ' ^ 


n+1 


'LOe 


'cJe 


Using Cauchy-Schwarz, Poincare and Young inequalities, we have 


f ~ ~n+l 

r Re -71 


n+1 


< L||j;||i2(^^)||7?. 11+2(0;,) 


n+1 


< rUp||j;||+2(^^)||V72. \\L‘^(eOe) 

< r(l — r)|| V7?. 


n+1, 


L-^{LOe) + 4(']_ _ ^^'-^Pll‘^e|lL2(o;,) 


ClWRe"^ 


As a matter of fact, r(l —r)||V'^ Ili 2 (t^^) can be absorbed by the heft-hand 

side of the energy estimate. The other source term can be treated in a very 
similar way. 

• linear terms 


J (jJe J 

Conventional arguments are the Holder inequality, the Sobolev injections like 
C with constant denoted C 5 , and Young’s inequality. For 

instance, the hrst term of Ce [TZ ), that is 7Z ■ Vmq, can be treated as 
follows: 


f . n+1 - n+1 

2r {n ■ Vuo) ■ n 
J cUe 


,, n+1,, o 

<2r\\n |Il4(^^)||Vmo||l2(^^) 

77 , I 1 

<2rCs\\S7'lZ ||i 2 (^^)||VMo||L 2 (o;,)- 


Under smallness assumptions on uq (and therefore on the data of the initial 
problem), the right-hand side of the above inequality can be absorbed by the 
left-hand side of the energy estimate. The other linear terms can be treated 
in a very similar way. 

• quadratic terms 


-Wee 


N 


QaiVlZ ,S 


n+1^ 


-n+l 
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Using the inequality 


Wee 


N 


— n -^n+l, -^n+l 

ga{VlZ ,s ): s 


n+l 


<Wee^||V7^1i2(,^)||5 
<Wee^C's||V^”|U2(,^)||V5 


^+1 II 


L'^iuJe) 


for e sufficiently small (namely Wee^ CsW'VlZ ||L 2 (tJs) < D), then the right- 
hand side of the above inequality can be absorbed by the left-hand side of 
the energy estimate. 

From the above considerations, we deduce the following estimate 
By Remark [9l we obtain 

^ , .v^ I 1 ^ , .v^ I 1 

IIhh-,) < 


( 10 ) 


|V7^ +1|5 "2 


Finally, in order to get the induction step on n, it is sufficient to guarantee that 

/V 9 -^n-^n 

Wee CsC < D to get the uniform estimate on [IZ , «S ). Note also that this 
condition is satisfied if £ is small enough. 

-^n-^n l/\ 

Step 2. We prove that (IZ , «S )„ is a Cauchy sequence in H [ue). Introducing 


n+1 ~ n+1 n 

iz .= iz -TZ , 


-^n+l ~n+l -^n 

s :=S -5 , 


QU+l _ QU+l _ QU^ 


we get by subtraction in (Run.) 

»r ~n n+1 ~ n, 

Re 7^ ■ V7^ + IZ ■ VIZ + Cr\lZ 


•'71+1 


n+1 


(1 - r)AlZ + VQ 
= div(5 


n+1 


n+1. 


n+1. 


and 


e^We[n VS +1Z VS +^„(V7^,5 )+ga{VlZ,S 

+ £f 

div'R^'^^ = 0. 


Performing an energy estimate consists, again, in multiplying the first equation 

^ n+1 -^n+1 

by 2r'TZ then integrate over cug, multiplying the second equation by S then 
integrate over and sum up the two contributions. We use the following estimates: 
• The first contributions is easily controlled, as 


£^Re 


— n n++ n+1 

iiz VIZ )-lZ 




= 0 . 


The second contribution satisfies: 


e^Re / {IZ VIZ ))-lZ 


n+1 


'CJe 
-.TV • 


■ n+1, 


< £ Re||7?. ||l 4 (^^)||V 7 ?. ||l 2 (^^)|| 7 ?. WL^^ico,) 

l”'+l in 


< r(l 


V7^ 


+ 


C £2A-2||y^ 


12 


' 4r(l—r)^ 

where we have used the estimate established in Eq. (Iiup and then Sobolev 
imbedding and Young’s inequality. 









VISCOELASTIC FLOWS IN A ROUGH CHANNEL: 


A MULTISCALE ANALYSIS 23 


• The third contribution satisfies 

J We 

where ||| ■ ||| denotes the operator norm from L‘^{us) to itself. Recalling 

the expression of it can be shown that |||£i"^^||| is arbitrarily small for 
sufficiently small data or e. 

• The other contributions can be treated in a straightforward way or by argu¬ 
ments similar to the previous ones. 

We thus deduce that, under smallness assumptions on the data and s, 

~ n+1 n ~n+l 

It means in particular that (7?. —IZ^S — ^ ) tends to 0 as n goes to +oc. 
Consequently, [IZ , «S ) is a Cauchy sequence in H (Ue). The sequence converges 
to some {IZ , «S ) in which, by means of construction, satisfies the following 

estimate: 

(11) IIVW + 115*111,.,-.) ^ 

Step 3. Letting n tend to -|-cxd we come to the conclusion that {IZ ,<S ) is the 
unique solution of the system satisfied by the lifted remainder. Therefore, the lifted 
remainder satisfies: 

IIV7?.||j;^ 2 -|- ||«S||j|fi Cs ^ 

which concludes the proof. □ 

Corollary 3. The asymptotic expansion is valid at any order. 


< 


14 "^^ I 




6. Algorithm and numerical procedure 


We show in this section that it is possible to effectively evaluate all the con¬ 
tributions of the ansatz. The only difficulty is to prove that the solution of each 
problem can be built using the only previous elementary solutions. For instance, let 
us prove that {uk,Pk, CTk), which is the solution of the problem pb^^\ only depends 
on {uj,pj,aj), j < k and on {Uj, Pj,Tij), j < k. This is not so obvious since the 
problem pb^^'^ calls for the use of a parameter a^, which seems to be related to a 
forthcoming problem PR(^+2) through the relationship 

Ok = lim [ F,^2{X,Y)dX. 

Y —>- + 00 

However we know prove that the definition of is a consistent: 

Proposition 6. Coefficients (see the definition o/pb*-^^ and satisfying 

Eq. ([H]) can he written as 


Oik 


lim / 

Y ->■- 1-00 Jr^d 


Re Mo • XxUk + (1 - r)A^Uk 

k—1 k 

^ ^ Ui * V xUj^—i ^ ^ Ui ’ (V xUk—i ~t“ V 

i=l i=l 
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Proof. The source contribution Fk +2 can be written as + ^k +2 with 

- Re Mo • V^f/fc + (1 - r)A^Uk + diVj,(Sfc) 

k-l 

^ ^ ' i.^^k+1—i T ^xUk—i) 

i=l 
k 

~ Re Ui ■ (Vt/fc+i-i + '^xUk-i + '\/uk-i) , 
i=l 

Fk +2 = - Rewo • Vt/fc+i - Rewfc • VUi + div(Sfc+i) 

+ 2(1 — r)'V X • Vt/fc+l ~ Va:-Pfc+ 1 . 

The first contribution -F ^2 using the already defined elementary solutions {{uj,pj, aj), 
j < k and {Uj, Pj, Sj), j < k), whereas the second ones Fk +2 using elementary solu¬ 
tions of problem and 

We now prove that lim / Fffr, = 0. For instance, we treat the first contribu- 
y^+oo Jja + 

tion (the four other contributions are similarly treated) 

d 

uo ■ Vt/fc +1 = - U^o^^'^dyUk+i. 

e=i 

Taking the X-average, using the periodicity we obtain 

f uo-VUk+i = ulf+^^dY( f Uk+l). 

J’fd V Jjd / 

Due to the behavior of the mean value ^k+i (see the Proposition [T]), we have 

lim / Mo • Vf/fc+i = 0 . 

y—)-+00 Jj'd 

To conclude the proof, it suffices to note that some contributions of F ^2 vanish 
too. □ 

The problem pb^^'> calls for the use of a parameter bk, which seem to be related to 
a forthcoming problem through the relationship 

bk= [ Xfe+i(X)dX-D-i / Gk+i{X,Y)dXdY. 

JT<i Jy<0 

However we know prove that the definition of bk is a consistent: 

Proposition 7. Coefficients bk (see the definition o/pb*-^^ and satisfying 

Eg. m only depend on the elementary solutions {uj,pj,aj)j^k and {Uj, Pj,'Ej)j<k- 

Proof. The source contribution Gk+i only depends on the already defined elementary 
solutions {{uj,pj, <Jj), j < k and (t/j, Pj, Sj), j < k). Only one term in the boundary 
contribution Nk+i does not depend on these elementary solution: = (VH ■ 

Vx)crfc_i|^-. Nevertheless, its X-average is equals to 0. □ 

The resulting application of Propositions | 6 ] and [7] leads to the following algorithm 
which allows us to build the approximation of the solution at any hxed order of 
precision: 
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Initialization: 

Compute ao (compatibility of problem at the top) 

Compute bo (compatibility of problem at the bottom) 

Compute (uo,po,ao) 

Compute (Ui, Fi, El) 

Iterative process on k: 

FOR fc = 1,+00, DO 

Compute Ofc (compatibility of problem PB*'^’*'^^ at the top) 
Compute hk (compatibility of problem PB*'^’*'^^ at the bottom) 
Compute {uk,Pk,(^k) 

Compute {Uk+i,Pk+i,Ek+i) 


END 


7. Concluding remarks 

7.1. Boundary conditions. The diffusive Oldroyd model is generally associated to 
Neumann boundary conditions for the elastic stress tensor [30], as suggested by the 
interpretation of the stress diffusion term as arising from the diffusion of polymeric 
dumbbells [18]. However some authors also considered Dirichlet boundary conditions 
|31) or mixed boundary conditions [2]. We emphasize that the method that has been 
developped in this paper readily adapts to the consideration of Dirichlet boundary 
conditions: the definition of the elementary problems has to be adapted and, in 
particular, the behaviour at inhnity of the boundary layer elastic tensor corrector 
is completely determined by an exponential decay towards a constant which can be 
counter-balanced by using a suitable elementary problem at next order. 

7.2. Influence of the curvature of the channel. The analysis of the roughness 
effects has been led in a particular geometry: a channel. When considering space- 
varying profiles of the boundary such as nozzles or more general converging-diverging 
profiles, additional coupling effects have to be taken into account. Although the 
methodology developped in this paper still applies, source terms have to be added in 
the elementary problems in order to include the scale effects at the macroscopic scale 
onto the boundary term which serves as a correction in the boundary layer problems. 
However, there is no additional difficulty from the mathematical point of view, 
although it tends to increase the complexity of the description of the elementary 
problems. 

7.3. The non-diffusive model. In many studies the standard Oldroyd model is 
considered without diffusion of the elastic stress tensor. The link between the stan¬ 
dard model and the diffusive model has been investigated from the numerical point 
of view by the authors [TB] by considering the vanishing diffusion process in the dif¬ 
fusive model. In the context of the roughness issue, the adaptation of our analysis is 
questionable when considering the standard Oldroyd model only. Although the two 
models are very close (at least formally in the regime D —)■ 0), the asymptotic expan¬ 
sion proposed in this paper does not apply to the case D = 0, even formally. This 
is due to the degeneracy of the boundary conditions: first, the loss of the boundary 
conditions in the standard model prevents us from developping the current strategy 
which is based upon the correction of the boundary terms at higher orders ; second. 
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ignoring the specific treatment of the boundary conditions that was done in the 

diffusive model provides inconsistent elementary problems. 
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